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Abstract. The nature of type la supernovae (SNe la) is still unclear. Employing Eggleton's stellar evolution code with the 
optically thick wind assumption, we systematically studied the He star donor channel of SNe la, in which a carbon-oxygen 
white dwarf accretes material from a He main-sequence star or a He subgiant to increase its mass to the Chandrasekhar 
mass. We mapped out the initial parameters for producing SNe la in the orbital period-secondary mass plane for various WD 
masses from this channel. According to a detailed binary population synthesis approach, we find that this channel can produce 
SNe la with short delay times (~100Myr) implied by recent observations. We obtained many properties of the surviving 
companions of this channel after SN explosion, which can be verified by future observations. We also find that the surviving 
companions from the SN explosion scenario have a high spatial velocity (>400 km/s), which could be an alternative origin 
for hypervelocity stars (HVSs), especially for HVSs such as US 708. 
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1. INTRODUCTION 

Type la supernovae (SNe la) play an important role in astrophysics, especially in the study of cosmic evolution. They 
have been applied successfully in determining cosmological parameters (e.g., Q.m and Ha; [1, 2]). They are also a key 
part of our understanding of galactic chemical evolution owing to the main contribution of iron to their host galaxies. It 
is generally believed that SNe la are thermonuclear explosions of carbon-oxygen white dwarfs (CO WDs) in binaries. 
However, there is still no agreement on the nature of their progenitors [3—6], and no SN la progenitor system before 
the explosion has been conclusively identified [7—9]. 

Over the past few decades, two families of SN la progenitor models have been proposed, i.e., the double-degenerate 
(DD) and single-degenerate (SD) models. Of these two models, the SD model is widely accepted at present. It is 
suggested that the DD model, which involves the merger of two CO WDs [10—12], likely leads to an accretion- 
induced collapse rather than to an SN la [13]. For the SD model, the companion could be a main-sequence (MS) star or 
a slightly evolved subgiant star (WD + MS channel), or a red-giant star (WD + RG channel) [14—26]. Observationally, 
there is increasing evidence indicating that at least some SNe la may come from the SD model [27—30]. 

Meanwhile, a CO WD may also accrete material from a He star to increase its mass to the Chandrasekhar (Ch) 
mass, which is known as the He star donor channel. Yoon and Langer [31] followed the evolution of a CO WD + He 
star system, in which the WD can increase its mass to the Ch mass by accreting material from the He star. It is believed 
that WD + He star systems are generally originated from intermediate mass binaries, which may explain SNe la with 
short delay times implied by recent observations [32—35]. 

In this contribution, we aim to investigate SN la birthrates and delay times of the He star donor channel and to 
explore the properties of the surviving companions after SN explosion. In Sect. 2, we describe the numerical code 
for the binary evolution calculations and the binary evolutionary results. We describe the binary population synthesis 
(BPS) method and results in Sect. 3. Finally, a discussion is given in Sect. 4. 



2. BINARY EVOLUTION CALCULATIONS 



We use Eggleton's stellar evolution code [36, 37, 38] to calculate the binary evolutions of WD + He star systems. 
The code has been updated with the latest input physics over the past three decades [39, 40, 41]. Roche lobe overflow 
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FIGURE 1. Regions in the orbital period-secondary mass plane for WD binaries (at the onset of the RLOF stage) that produce 
SNe la for initial WD masses of 0.865,0.9,1.0,1.1 and 1.2 M . (From Wang et al. [45]) 



(RLOF) is treated within the code described by Han et al. [42]. We set the ratio of mixing length to local pressure 
scale height, a = l/H p , to be 2.0. The opacity tables are compiled by Chen and Tout [43]. In our calculations, He star 
models are composed by He abundance Y = 0.98 and metallicity Z = 0.02. 

Instead of solving stellar structure equations of a WD, we use the optically thick wind model [14] and adopt the 
prescription of Kato & Hachisu [44] for the mass accumulation efficiency of He-shell flashes onto the WD. We have 
calculated about 2600 WD + He star systems, and obtain a large, dense model grid (for details see [45]). In Fig. 1, we 
show the contours at the onset of RLOF for producing SNe la in the logP 1 — M\ plane for various WD masses (i.e., 
M^, D = 0.865, 0.90, 1.0, 1.1 and 1.2M©), where P 1 and M\ are the orbital period and the mass of the He companion 
star at the onset of RLOF, respectively. The region almost vanishes for M^, D =O.865M , which is then assumed to 
be the minimum WD mass for producing SNe la from this channel. If the parameters of a WD + He star system at 
the onset of the RLOF are located in the contours, an SN la is then assumed to be produced. Thus, these contours 
can be expediently used in BPS studies to investigate the He star donor channel of SNe la. The data points and the 
interpolation FORTRAN code for these contours can be supplied on request by contacting BW. 



3. BINARY POPULATION SYNTHESIS 

To obtain SN la birthrates and delay times of the He star donor channel, we performed a series of Monte Carlo 
simulations in the BPS study. In the simulation, by using the Hurley's rapid binary evolution code [46, 47], we 
followed the evolution of 4 x 10 7 sample binaries from the star formation to the formation of the WD + He star 
systems according to the SN la production regions (Fig. 1) and three evolutionary channels (i.e., the He star channel, 
the EAGB channel, and the TPAGB channel; for details see [48]). Here, we adopt the standard energy equations to 
calculate the output of the common envelop (CE) phase (see [48]). 

In the BPS study, the primordial binary samples are generated in the Monte Carlo way and a circular orbit is assumed 
for all binaries. We adopt the following input for the simulation (for details see [48]). (1) The initial mass function 
of Miller and Scalo [49] is adopted. (2) The mass-ratio distribution is taken to be constant. (3) The distribution of 
separations is taken to be constant in log a for wide binaries, where a is the orbital separation. (4) We simply assume 
a constant star formation rate (SFR) over the past 15 Gyr or, alternatively, as a delta function, i.e., a single starburst. 

In Fig. 2, we show the evolution of SN la birthrates with time for the He star donor channel. The left panel represents 
Galactic birthrates of SNe la by adopting Z = 0.02 and SFR = 5M©yr _1 . The simulations give Galactic SN la birthrates 
of ~0.3 x 10~ 3 yr _1 . The right panel displays the evolution of SN la birthrates for a single starburst with a total mass 
of 10 n M R . In this panel, we see that SN la explosions occur between ~45Myr and ~ 140 Myr after the starburst, 
which could explain SNe la with short delay times implied by recent observations [32—35]. 

The companion in the SD model would survive in the SN explosion and potentially be identifiable [50, 51]. We 



a OT >>- 0.5 



o 




0^=0.5 
a„X=1.5 



log(«/yr) 



log(t/yr) 



FIGURE 2. The evolution of SN la birthrates with time for the He star donor channel. Left panel: the evolution of Galactic SN 
la birthrates for a constant Pop I SFR (5M B yr -1 ), where the solid and dashed curves show the results of different CE ejection 
parameters with a ce A = 0.5 (solid) and a ce A = 1.5 (dashed), respectively. Right panel: similar to the left panel, but for a single 
starburst with a total mass of 10 11 M Q . (From Wang et al. [48]) 
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FIGURE 3. Left panel: the distribution of properties of companions in the plane of (V or b, M| N ) at the current epoch, where V or b 
is the orbital velocity and m| n the mass at the moment of SN explosion. Right panel: similar to the left panel, but in the plane of 
(logT e ff, logg), where T e ff is the effective temperature of companions and logg the surface gravity. (From Wang & Han [7]) 
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FIGURE 4. Left panel: similar to Fig. 3, but in the plane of (logP , Af| ), where P is the orbital period. Right panel: similar 



to Fig. 3, but in the plane of (V rot , M| N ), where V rot is the equatorial rotational velocity of companions. (From Wang & Han [7]) 
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FIGURE 5. The distribution of the spatial velocity with a ce A = 0.5. (From Wang & Han [7]) 



obtained the distributions of many properties of the surviving companions of this channel at the moment of SN 
explosion. Figs 3 and 4 shows the distributions of the masses, the orbital velocities, the effective temperatures, the 
surface gravities, the orbital period and the equatorial rotational velocity of companions at the moment of SN explosion 
(for the detailed discussion see [7]). These distributions may be helpful for identifying the surviving companions of 
SNe la. The simulation also shows the initial parameters of the primordial binaries and the WD + He star systems that 
lead to SNe la, which may help to search for potential SN la progenitors. 

In Fig. 5, we show the current epoch distribution of the spatial velocity for the surviving companions from this 
channel (for details see [7]). We see that the surviving companions have high spatial velocities (>400km/s), which 
almost exceed the gravitational pull of the Galaxy nearby the Sun. Thus, the surviving companions from the SN 
explosion scenario could be an alternative origin for hypervelocity stars (HVSs), which are stars with a velocity so 
great that they are able to escape the gravitational pull of the Galaxy. 



4. DISCUSSION 

The simulations give Galactic SN la birthrate of ~0.3 x 1CT 3 yr _1 for the He star donor channel, which is lower than 
that inferred observationally (i.e., 3 — 4 x 10~ 3 yr _1 ; Cappellaro and Turatto [52]). This implies that the He star donor 
channel is only a subclass of SN la production, and there may be some other channels or mechanisms also contributing 
to SNe la, e.g., WD + MS channel, WD + RG channel or double-degenerate channel (see [53]). 

Hachisu et al. [54] recently studied the mass-stripping effect on a MS or slightly evolved companion star by winds 
from a mass-accreting WD. The model can also provide a possible ways of producing young SNe la, but the model 
depends on the efficiency of the mass-stripping effect. We also find that the model produces very few young SNe 
la according to a detailed BPS approach. Thus, we consider the He star donor channel as a main contribution to the 
formation of young SNe la. 

Based on the optically thick wind assumption, Wang & Han [55] recently calculated about 10 000 WD + He star 
systems and obtained SN la production regions of the He star donor channel with different metallicities. For a constant 
star-formation galaxy, they found that SN la birthrates increase with metallicity. If a single starburst is assumed, SNe 
la occur systemically earlier and the peak value of the birthrate is larger for a high Z. We also note that Liu et al. [56] 
recently investigated the effects of rapid differential rotation on the accreting WD, and found that the highest mass of 
the accreting WD at the moment of SN la explosion is 1.81M© for the He star donor channel, which may provide a 
way for the formation of super-Ch mass SNe la in the observations. 

Massive WD + He star systems are candidates of SN la progenitors. HD 49798/RX J0648. 0-4418 is an evidence 
of the existence of massive WD + He star systems. Based on the data from the XMM-Newton satellite, Mereghetti 
et al. [57] recently derived the masses of the two components. The corresponding masses are 1.50±0.05M© for HD 




FIGURE 6. Results of binary evolution calculations with initial masses of two components and the orbital period as well similar 
to the binary system HD49798/RX J0648. 9-4418. In the left panel, the solid, dashed and dash-dotted curves show the mass-transfer 
rate, the mass-growth rate of the WD and the mass of the WD varying with time after HD 49798 fills its Roche lobe, respectively. 
In the right panel, the evolutionary track of the He donor star is shown as a solid curve, and the evolution of the orbital period is 
shown as a dash-dotted curve. The solid star represents the current position of HD 49798. Dotted vertical lines in these two panels 
and open stars in the right panel indicate the position where the WD is expected to explode as an SN la. (From Wang & Han [58]) 



49798 and 1.28±O.O5M for the WD. According to a detailed binary evolution model, Wang & Han [58] found that 
the massive WD can increase its mass to the Ch mass in future evolution (see Fig. 6). Thus, HD 49798 with its WD 
companion is a likely candidate of SN la progenitors. V445 Pup is another candidate of massive WD + He star system, 
in which the mass of the WD is more than 1.35M©, and the mass of the He star is more than 0.85M© [59, 60]. 
However, we still do not know the orbital period of the binary system and the mass of the He donor star so far. This 
needs further observations of V445 Pup after the dense dust shell disappears. 

US 708 is an extremely He-rich sdO star in the Galaxy halo, with a heliocentric radial velocity of +708 ± 15km/s 
(Hirsch et al. [61]). We note that the local velocity relative to the Galatic center may lead to a higher observation 
velocity for the surviving companions, but this may also lead to a lower observation velocity. Considering the local 
velocity nearby the Sun (^220km/s), we find that ^30% of the surviving companions may be observed to have 
velocity V > 700km/s for a given SN ejecta velocity 13500 km/s. In addition, the asymmetric explosion of SNe la 
may also enhance the velocity of the surviving companions. Thus, a surviving companion in the He star donor channel 
may have a high velocity like US 708. In future investigations, we will employ the Large sky Area Multi-Object fiber 
Spectral Telescope (LAMOST) to search the HVSs originating from the surviving companions of SNe la. 
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